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Abstract — In this paper, we consider the demodulation and 
equahzation problem of differential Impulse Radio (IR) Ultra- 
WideBand (UWB) Systems with Inter-Symbol-Interference (ISI). 
The differential IR UWB systems have been extensively discussed 
recently [1], [2], [3], [4], [5], [6], [7]. The advantage of differential 
' IR UWB systems include simple receiver frontend structure. One 
challenge in the demodulation and equalization of such systems 
with ISI is that the systems have a rather complex model. The 
input and output signals of the systems follow a second-order 
Volterra model [7]. Furthermore, the noise at the output is data 
dependent. In this paper, we propose a reduced-complexity joint 
demodulation and equalization algorithm. The algorithm is based 
on reformulating the nearest neighborhood decoding problem 
into a mixed quadratic programming and utilizing a semi- 
definite relaxation. The numerical results show that the proposed 
demodulation and equalization algorithm has low computational 
complexity, and at the same time, has almost the same error 
probability performance compared with the maximal likelihood 
decoding algorithm. 

I I. Introduction 

Ultra- WideBand (UWB) communication systems have at- 
tracted much attention recently. The UWB communications 
have the advantages of robustness due to multi-path diversity, 

' low possibilities of intercept and high location estimation 
accuracy. UWB systems are favorable choices for short range 
high bit rate communications or medium-to-long range low bit 

\ rate communications. For example, UWB systems have been 
considered for video communications in Wireless Personal 
Area Networks (WPAN). In this case, the transmission rates 
can be as high as 400M bits per second. UWB communica- 
tion systems have also been considered for Wireless Sensor 
Networks (WSN) as a low-power and low-cost solution. The 
FCC (US Federal Communications Commission) has recently 
approved the use of UWB communications and allocated a 
spectrum range of 7.5 GHz for UWB communications. 

A communication system is considered to be a UWB 
system, if the system's bandwidth spans more than 1.5 GHz, or 
25% of the center frequency. The UWB systems transmit data 
by sending pulses, each with very small time duration. For one 
transmitted pulse, a large number of replicas of the same pulse 
are received at the receiver side due to multi-path. The number 
of resolvable multi-paths can be as high as more than 100 
as shown in [8]. As a consequence, multi-path diversities are 
automatically achieved. However, accurate channel estimation 
can be quite complex and difficult. 



The existing approaches for UWB communications include, 
Direct-Sequence (DS) UWB, Multi-Band (MB) UWB, and 
low-complexity non-coherent Impulse Radio (IR) UWB sys- 
tems. The DS-UWB systems use direct sequence spreading 
technique to convert information signals into wide-band sig- 
nals, [9], [10]. Under the condition that the channel estimation 
is accurate, the RAKE receiver is the optimal demodulation 
scheme. However, the channel estimation for UWB channels 
is difficult and complex. Without the information about the 
correct RAKE weights, the systems suffer a performance loss 
by using sub-optimal RAKE structures (for example, equal 
weight combining). 

MB-UWB systems are recently proposed and discussed in 
[11], [12], [13]. The MB-UWB systems use the Orthogonal 
Frequency Division Multiplexing (OFDM) technology. The 
advantage of MB-UWB systems include higher achievable 
bit rates, flexibility in spectrum occupation, good coexistence 
with narrow band communications. The disadvantages include 
complex architectures, and high power consumption. The third 
class of UWB systems is the non-coherent IR UWB systems. 
In such systems, complete channel estimation is not required. 
Therefore, the channel estimation constraint is greatly relaxed. 

In this paper, we consider a low-complexity non-coherent 
IR UWB system - the differential IR UWB system proposed 
in [4]. In the differential IR UWB systems, the transmitted 
information is differentially encoded. At the receiver side, 
a low-complexity Autocorrelation (AcR) receiver is adopted. 
The decoding decision variables are autocorrelations 

f \{t)r{t + S)dt, (1) 

J to 

where, r{t) is the received signal, and S is the time difference 
between two consecutive pulses. The integral can be imple- 
mented either in the analog domain or in the digital domain. 
In both cases, the decoder architecture is largely simplified. 

One problem of the AcR receiver is that the transmitted 
messages and the receiver decoding decision variables follow 
a nonlinear second-order Volterra model, especially when 
Inter-Symbol-Interference (ISI) is present in the systems [7]. 
The maximal-likelihood sequential decoders can be adopted, 
however their computational complexities generally grow ex- 
ponentially with the length of delay spread. 

In this paper, we propose a reduced-complexity demod- 
ulation and equahzation algorithm. The algorithm is based 
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on a reformulation of the nearest neighborhood decoding 
problem into a mixed quadratic programming and a Semi- 
Definite Programming (SDP) relaxation. The computational 
complexity of the proposed algorithm grows only polynomi- 
ally with respect to the block length and is independent of 
the length of delay spread. We show by simulation results 
that the performance loss caused by the proposed sub-optimal 
demodulation algorithm is negligible. 

SDP relaxation has been previously adopted to solve de- 
coding problems and combinatorial optimization problems. In 
[14], an approximation algorithm for maximum cut problem 
based on SDP relaxation has been proposed. Detection algo- 
rithms for MIMO channels based on SDP relaxation have also 
been proposed in [15], [16], [17], [18], [19]. For interested 
readers, a review of SDP optimization can be found in [20]. 

The rest of this paper is organized as follows. In Section 
im we describe the system model. We present the proposed 
demodulation and equalization algorithm in Section Nu- 
merical results are presented in Section |IV] Conclusions are 
presented in Section |V] 

Notation: We use the symbol S to denote the set of 
symmetric matrices. Matrices are denoted by upper bold face 
letters and column vectors are denoted by lower bold face 
letters. We use A ^ to denote that the matrix A is positive 
semi-definite. The symbol ® is used to denote the Kronecker 
product. We use Aij to denote the element of the matrix A 
at the i-th row and j-th column. We use Oi to denote the i-th 
element of the vector a. We use A^ and to denote the 
transpose of the matrix A and the vector a respectively. We 
use tr{A) to denote the trace of the matrix A. The function 
sign(-) is defined as. 



sign(x) 



1, if a; > 
— 1 , otherwise 



(2) 



II. System Model 

We assume that the message is transmitted in a block by 
block fashion. The transmitted signal in one block is 

Wb-l Np-l 

s(t)= J2 J2 a.,[n]iu {t - U[n]) (3) 

n=0 i=0 

where w{t) is the transmitted pulse, ai [n] is the pulse polarity 
for the i-th pulse of the n-th symbol, ti [n] is the pulse time for 
the i-th pulse of the n-th symbol. Each block has Nf, symbols, 
and each symbol corresponds to Np pulses. 

Denote the data symbol by d[n] e { — The data 
symbols are differentially encoded as. 



ai\n 



aArp_i[n - - l]6jVp-i, if i = 

ai-i[n]d[n]bi^i, otherwise 



(4) 



where, bo,bi, . . . ,bN is the pseudo-random amplitude code 
sequence, bi G { — The pulse time 



U [n] = nTs + Ci 



(5) 



where Tg is the symbol duration, Cj is the relative pulse timing. 
The relative pulse timeing Ci is related to the pseudo-random 



delay hopping code {D,}, 

Ts+co~CNp-i, if Np-l 



D, 



Ci+l - Ci, 



otherwise 



(6) 



The pseudo-random amplitude code and delay hopping code 
are used to facilitate multiple access. 
The received signal is, 

^W=£ 51 "»["]ff(i-i»W)+"W, (7) 

where, g{t) is the channel response for the pulse w{t), n{t) is 
the noise. The receiver front end is shown in Fig. [T] Denote 
the decoding decision variable for the n-th symbol by z[n\, 

<.t,[n]+T/ 

y,[n] = / r{t)r{t + D,)dt, (8) 

Jtiln] 



(9) 



where, T/ is the integral time. 

r(t) 
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Fig. 1. Block diagram of the autocorrelator receiver 

Let us define 

Ig{h,t2;T) = f g{t)g{t + T)dt 
Jti 

Denote the data vector by d, 

d ~ [do, di, . . . , C?Ar^_i]"^. 

Define the column vector, 

a = [ao[0],ai[0], . . . , ao[n], ai[n], . . .,aNp^i[Nb - 1]] 

Neglecting noise, we have 

yi[n\ = a^Ai[n\a. 



(10) 
(11) 

T 

(12) 
(13) 



In the above equation, Ai[n] is a matrix, such that the {n'Np + 
i' + 1, n"Np + i" + 1) element is 

Ig {U[n] ~ te[n']Mn\ - U'[n'] + Ti-U>[n'] - U»[n"] + D,) . 

(14) 
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Finally, the decoding decision variables can be written as, 

Np-l 

z[n] — ^ bia'^Ai[n]a = a^B[n]a, (15) 



1=0 



where, B[n\ = J2^=o ^ i'i-^ii'n]- 



(16) 



The vector a can be written as, 

a = Q{r + Pd). 
In the above equation, Q is a diagonal matrix, 

Q=diag[l,5o,5o6i,...], (17) 



The matrix P, 



fc-2 

[Qhk = n ^modATp- 



P = /at, S, 



(18) 



(19) 



where, / is an identical matrix, and s is a vector with length 
Np of alternating 0, 1, 



The vector r is. 



s= [0,1,0,1,...,!]^ 



(20) 



(21) 



where i and i n are the all one column vectors with length 
Nb and Np respectively. 

With the above notation, the second-order Volterra model 
of the system is, 

z[n] = {r + Pd)^Q^B[n]Q{r + Pd) + noise terms, (22) 

where the noise terms are data dependent as shown in [7]. 



III. Joint Demodulation and Equalization 
Algorithm 

In this section, we present the proposed demodulation and 
equalization algorithm. The algorithm is obtained by formu- 
lating the demodulation problem as a nearest neighborhood 
decoding problem, reformulating into mixed quadratic pro- 
gramming, and using SDP relaxation. 

In the first step, we formulate the demodulation problem as 
a nearest neighborhood decoding problem as follows. 

min {z[n] - {r + PdfQ^B[n]Q{r + Pd)}^ (23) 



subject to dn G { — 1,1}. 



(24) 



Note that the nearest neighborhood decoding is not the max- 
imal likelihood decoding in the considered scenario, because 
noise is signal dependent. 



The above nearest neighborhood decoding problem can be 
reformulated as a mixed quadratic programming by introduc- 
ing auxiliary variables s„, n = 0, . . . , N}, — 1. 



n=0 



(25) 



(26) 



subject to 

Sn = z[n]-{r + PdfQ^B[n]Q{r + Pd), (27) 



dn e {-1,1}. 



(28) 



Now, we claim that the above mixed quadratic programming 
is equivalent to the following matrix optimization problem. 



Nt + l 

min ^ Ur, 



(29) 
(30) 



n=2 

subject to 
E/i „ = z[n - 2] - r^Q^B[n - 2]Qr 

-r^Q^B[n 2]QPd' 

-r^Q^B[n-2fQPd' (31) 

- tr {D'P^Q^B[n - 2]QP} , 

for n = 2, . . . , iVb + 1 
Ui,i = 1, (32) 
Un,n = 1, for n = iVf, + 2, . . . , 2iV6 + 1, (33) 

d' = [Ui,N,+2,.-.,Ui,2N, + lf, (34) 

U eS, (35) 
U hO, (36) 
U has rank one, (37) 

where, U denotes a matrix of size 2Nt, + 1 by 27Vb + 1, 
D' denote the sub-matrix of U formed by selecting the last 
A^f, rows and columns, and the optimization variables are the 
elements of the matrix U. Because the matrix U has rank one, 
is symmetric and positive semi-definite, it is well known [20] 
that there exists a vector u, such that 



U = uu^. 



(38) 



If we further assume that Ui = 1, then the vector u is unique. 
In addition, there is an one-to-one correspondence between the 
solution of the mixed quadratic programming and the solution 
of the matrix optimization problem. 



(39) 



Therefore, the mixed quadratic programming is equivalent to 
the matrix optimization problem. 

In the matrix optimization problem, all objective function 
and constraints are convex except the rank one constraint in 
Eq. [37] If the rank one constraint is relaxed, then we obtain a 
convex relaxation (SDP relaxation). The convex optimization 
problem can then be efficiently solved by polynomial-time 
algorithms and softwares, for example, by using the SeDuMi 
package [21]. Previous research has shown that such SDP 
relaxations are tight approximations to the original problems 
[14]. Near-optimal solutions of the original problems can be 
obtained from SDP relaxations by random rounding. 
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Finally, the proposed joint demodulation and equalization 
algorithm consists of two steps. In the first step, the following 
SDP relaxation problem is solved. 

Nb + l 

min ^ Un,n (40) 

n=2 

subject to (41) 
Z7i „ = z[n - 2] - r^Q^B[n - 2]Qr 
-r^Q^B[n-2]QPd' 

- r^Q^B[n - 2fQPd' (42) 

- tr {D'P'^Q'^B[n - 2\QP) , 
for n = 2, . . . , TVf, + 1 

Ui,i = 1, (43) 
Un,n = 1, for n = iVb + 2, . . . , 2iVb + 1, (44) 
= [t^i,JV,+2,...,C/i,2iv.+i]^, (45) 
U eS, (46) 
U>0, (47) 
D' is the submatrix of U formed by selecting the last Ni, 
rows and columns. (48) 

In the second step, the demodulation decision is made by 
thresholding. 



dn = sign(i7i,„+jVt,+2)- 



(49) 



IV. 



Numerical Results 

In this section, we present simulation results for the pro- 
posed demodulation and equalization scheme. We assume that 
the channel can be modeled by the S-V model [22]. The 
received signal for each transmitted pulse w{t) is. 



9(t) 



E 

i=i 



ajw{t — Sj). 



(50) 



In the above equation, w{t) is the second derivative Gaussian 
monocycle. 



w{t) = [l-47r(t/r™)'] exp|-27r(t/r„)'} 



(51) 



where t„i = 0.2877 nanosecond. N„i is the total number of 
multiple paths, aj and 6j are amplitude and delay of the j-th 
path. 

We assume that the delays of the paths follow the Poisson 
process with the expected interval between two consecutive 
paths being 10 nanoseconds. The amplitude is Raleigh dis- 
tributed, such that the expectation of the amplitude aj is 
exp(— (Jj/Te), where Tg = 20 nanoseconds. We assume that 
the amplitudes and delays aj, 5j vary slowly, so that the 
matrices B[n] can be accurately estimated (for example, by 
using pilot signals), and considered perfectly known at the 
demodulator 

For the transmitted signal, we assume that each block 
has Nb = 10 symbols and each symbol corresponds to 
Np = 4 pulses. The pseudo-random delay hopping code is 
[1.7, 1.9,2.1,2.3]. The symbol duration Tg — 8 nanoseconds. 
The integral time Tj ^ Tg = 8 nanoseconds. 



We illustrate the bit error probability of the proposed 
demodulation and equalization algorithm in three different 
cases. In the first case, we assume that the delay spread extends 
over a range of 200 nanoseconds. Therefore, there exists severe 
non-linearity in the system. The bit error probability of the 
proposed scheme for this case is shown in Fig. |2] The bit 
error probabihty of the maximal likelihood detection algorithm 
is also plotted. 

In the second case, we assume that the delay spread extends 
over a range of 30 nanoseconds. The non-linearity is mild in 
this case. The bit error probabilities of the proposed scheme 
and the maximal likelihood decoding are shown in Fig. [3] In 
the third case, we assume that there is no ISI. And the bit 
error probabilities are shown in Fig. |4] From all these results, 
we conclude that even though the proposed scheme is sub- 
optimal, the performance loss is negligible. 



10" 



10" 



Bit error probability 




maximal likelihood decoding 

convex optimization based 



14 16 18 20 22 

energy per bit to noise power spectral density ratio in db 



Fig. 2. Bit eiTor probabilities in the case of severe non-linearity. The length 
of delay spread is 200 nanoseconds. The solid curve represents the bit error 
probabilities of the maximal likelihood detection algorithm. The dashed curve 
represents the bit error probabilities of the proposed reduced complexity 
detection algorithm. The X-axis shows energy per bit to noise power spectral 
density ratio Et /No in dB. 



V. Conclusion 

In this paper, we propose a convex optimization based 
demodulation and equalization algorithm with low complexity 
for the differential IR UWB systems. The complexity of the 
proposed algorithm grows polynomially with respect to the 
blocklengths, and is independent of the length of delay spread. 
Even though the proposed algorithm is sub-optimal, we show 
by simulation results that the performance loss is negligible. 
The proposed demodulation and equalization algorithm is a 
near-optimal algorithm with significantly reduced computa- 
tional complexity. 
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